Symmetry breaking of the zero energy Landau level in bilayer graphene 
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The quantum Hall effect near the charge neutrality point in bilayer graphene is investigated in high 
magnetic fields of up to 35 T using electronic transport measurements. In the high field regime, the 
eight-fold degeneracy in the zero energy Landau level is completely lifted, exhibiting new quantum 
Hall states corresponding filling factors v =0, 1, 2, & 3. Measurements of the activation energy gap 
in tilted magnetic fields suggest that the Landau level splitting at the newly formed v =1, 2, & 3 
filling factors are independent of spin, consistent with the formation of a quantum Hall ferromagnet. 
In addition, measurements taken at the v — charge neutral point show that, similar to single layer 
graphene, the bilayer becomes insulating at high fields. 

PACS numbers: 73.63.-b, 65.80.+n, 73.22.-f 
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The unique chiral nature of the carrier dynamics in 
graphene results in a novel integer quantum Hall (QH) 
effect that distinguishes these atomically thin graphitic 
materials from conventional 2-dimensional systems [l|, 0, 
3]. The chiral nature of carriers in single layer graphene 
(SLG) and bilayer graphene (BLG) result in unevenly 
spaced Landau levels (LL) including a distinctive level 
located precisely at the particle-hole degenerate zero en- 
ergy. While the energy spacing of SLG LLs scales as the 
square root of the LL index n and the square root of the 
field B, BLG levels have energies given by 0] 



±ftw\/ n(n — 1), 
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where lu = eB/m is the cyclotron frequency with the 
BLG band mass m ~ 0.05 m e , and n is a non-negative 
integer representing Landau orbit index in each layer. 
While the valley and spin degrees of freedom lead to a 
4-fold degeneracy in graphene LLs, the additional orbital 
degeneracy for the n =0 & 1 indices in BLG result in a 
8-fold degenerate, zero energy level that is unprecedented 
in LL physics dHH. 

At high fields, the decreasing radius of the cyclotron or- 
bits gives rise to increasing electron-electron interactions 
which can perturb the degenerate LLs. Experiments on 
SLG have demonstrated this high-field symmetry break- 
ing through the appearance of new QH states at the fill- 
ing factor sequences v = 0,±1,±4 jy, 0, S|- The precise 
nature of the field-dependent mechanisms that lift the 
degeneracies is still under experimental and theoretical 
debate 0, [lfj. In particular, it has recently been ob- 
served that in SLG the v = QH state becomes increas- 
ingly insulating at higher fields 11, 12j]. Similar to SLG, 
the enhanced interactions under high magnetic field is 
expected to lift the 8-fold zero energy LL degeneracy in 
BLG. A number of theoretical predictions involving un- 
usual collective excitations have been proposed to occur 
in this particle-hole symmetric LL as its degeneracy is 
broken jl, H Q 

Il5j . but the experimental observation 



of this broken symmetry is yet to be observed. 

In this Letter we present transport measurements 
demonstrating the hierarchy of the splittings for the zero- 
energy LL in bilayer graphene as the external magnetic 
field is increased. Activation energy measurements of 
these splittings demonstrate energy spacing smaller than 
the bare electron-electron interaction energy and Zccman 
energy. Tilted field measurements indicate that the de- 
generacy breaking at filling factors above the charge neu- 
tral point are independent of spin, consistent with the 
formation of a quantum Hall ferromagnet. While a quan- 
tum Hall plateau in the transverse resistance is observed 
at the charge neutral point, the longitudinal resistance at 
this point is seen to diverge to an insulating state with 
increasing field. 

The graphene samples used in this work are deposited 
on SiO2(300 nm)/Si substrates by mechanical exfoliation 
techniques from bulk single crystals [l6| . The number of 
graphene layers are identified by optical contrast, cross- 
checked using Raman spectroscopy fl7j j and measure- 
ments of the high-field LL spectrum. Au/Cr electrodes 
for the graphene are fabricated by conventional electron 
beam lithography followed by metal evaporation. A gate 
voltage V g is applied to the degenerately doped Si sub- 
strate to control the carrier density n c according to the 
relation n c = C g (V g — Vh)/e, where the areal gate ca- 
pacitance is C g = 7.1 x 10 10 e/V • cm 2 and Vd is the gate 
voltage corresponding the charge neutrality point. The 
lower inset of Fig. la shows an optical microscope image 
of a typical BLG device used in this study. This device 
has a mobility as high as ~ 1 x 10 4 cm 2 /V • s measured 
at the carrier density n/j = 4 x 10 12 cm~ 2 . Since the elec- 
trodes of this device are configured in non-ideal Hall bar 
geometry, the magnetoresistance R xx and Hall resistance 
R xy were obtained following the van der Pauw method 
with symmetric (R xx ) and anti-symmetric (R xy ) averag- 
ing over data from positive and negative magnetic fields. 
The Hall conductivity a xy is then computed from R xx 
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FIG. 1: (a) Hall conductivity a xy , as a function of gate voltage 
V g at T=1.4K at different magnetic fields: 9, 15, 20, 25, 30, 
and 35 T. Upper inset: R xx (in black) and R xy (in red) as V g 
varies at B=35 T. Lower inset: optical microscope image of 
a BLG device used in this experiment, (b) Schematic of the 
zero energy LL hierarchy in bilayer graphene at high magnetic 
field. 

and R xy . 

Fig. 1(a) shows R xx and R xy (upper inset) and a xy . 
As the back gate voltage alters the carrier density, QH 
plateaus in R xy with corresponding zeros in R xx are 
observed. Consequently, the calculated <r xy shows well 
quantized plateau values at ~-tt with the integer filling 
factor v. In the low magnetic field regime [B < 10 T), 
QH states corresponding to v = ±4, ±8 are present, as 
were previously seen in Ref. As B increases, how- 
ever, new QH states emerge, as evidenced by additional 
QH plateaus at v = & 2 and then v = 1 & 3 at higher 
fields. 

In the high magnetic field regime, B > 20 T, the new 
set of QH states demonstrate the full 8-fold degeneracy 
lifting in the zero energy LL. On close inspection of a xy 
in different magnetic fields, we can construct the hier- 
archical evolution of these new QH states. As the field 
increases from the low magnetic field regime, the v = 
and 2 states start to develop at 15 T and are fully evolved 
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V g /(B tot cos(9)) 

FIG. 2: (a) Magnetoresistance R xx as a function of normal- 
ized back gate voltage V g / (Btotcos(8)) around v = 2 and 
v = 3, at four different total magnetic fields with the same 
perpendicular field. Upper inset: Schematic diagram of tilted 
field; Lower inset: R xx vs. V g / (Btotcos(6)) at v = 1 state 
under the same condition as in the main figure. (b)R xx as a 
function of normalized gate voltage with a zero tilting angle 
at different total fields, with the same scale of Fig2(a). Lower 
inset: R xx vs. V g / (B to tcos(8)) at u = 1 state. All data are 
taken at 1.5 K. 



by 20 T; while the v = 1 and v = 3 states start to develop 
at 20 T and are fully evolved by 25 T. This hierarchical 
appearance of the new QH states is in accordance with 
the sequential symmetry breaking of the zero energy LL 
degeneracy as depicted in Fig. 1(b) and suggests that 
different symmetry-breaking process are relevant as B 
increases. 

Mechanisms for breaking the LL degeneracy include 
disorder, lattice strain, and charged impurities. The 
magnitude of the last of these can be estimated by assum- 
ing that the offset of the Dirac point gate voltage Vd from 
zero energy reflects a bilayer capacitively charged across 
the d = 0.34nm layer spacing. From this model our mea- 
sured I Vd\ ~ 2 V corresponds to a maximum charging en- 
ergy of ~ 2 meV, small than, but of the same order as the 
high-field (B > 20T) bare Zeeman energy. In addition to 
the above mentioned graphene impurities, there are two 
field-dependent factors that can lead to the lifting of the 
eight degeneracies of this LL: Zeeman splitting of spin 
in magnetic field and electron-electron interactions that 
increase with field as the radius of the cyclotron orbits de- 
crease. Zeeman splitting is given by A z = g/isB, where 
fiB is the Bohr magneton, and g is the gyromagnetic fac- 
tor for the carriers in BLG. The Coulomb interactions 
between electrons are given by e 2 /els, where e ~ 4 is 
the dielectric constant for graphene and Ib — \/h/eB 
is the magnetic length. In order to produce full 8-fold 
symmetry breaking, these mechanisms must lift not only 
the spin and valley degeneracies, but the n = & 1 LL 
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orbital degeneracy. Barlas et al. have suggested that the 
exchange term from Coulomb interactions between elec- 
trons in these degenerate orbital states leads to a QH fer- 
romagnetic state where the 8-fold symmetry is lifted by 
this exchange interaction. A Hund rule-like hierarchical 
symmetry breaking in the zero energy LL is predicted 
with this exchange-enhanced Zeeman splitting followed 
by the spin-independent valley and orbital splitting [5]. 

In order to test out this hierarchical degenerate lifting 
and the role of electron-electron interactions, we first note 
that sublattice or interlayer splittings that are associated 
solely with the Coulomb interactions between electrons 
localized in LL orbits should depend only on the nor- 
mal magnetic field (B±) but not the in-plane field (B\\). 
Whereas spin-splitting includes a Zeeman component de- 
termined by the total field B to t - 



B%. Experi- 
mentally, this idea can be tested out by measuring the 
samples in a series of tilted fields where we can exam- 
ine the R xx minima in different B± and B\\ by tuning 
the tilting angle 9 and B tot . Fig. 2 displays the change 
of R xx as a function of the gate voltage normalized by 
B± = B to t cos 9, in two different experimental conditions. 

In Fig. 2(a), we first fix B p — 20 T and then vary 
Bh. For v — 2,3 QH states, -R™", the minima of R xx 
corresponding to the QH zero magnetoresistance, display 
less than ~ 10% of variation, indicating no significant 
changes are induced by applying B^ . For comparison, in 
Fig. 2(b), we show R xx as a function of B± with fixed 
9 = 0, that is By = and B ± = B tot . In this case R™ n 
for v = 2 and 3 increases by ~ 40% as B± decreases 
from 31 T to 20 T, a factor of 4 larger change than the 
change of Fig. 2(a) where B to t decreases from 31 T to 
20 T while B ± was kept to 20 T. The fact that E™ n 
for v = 2 and 3 shows little to no dependence on the in- 
plane field strongly suggests that both these QH states 
are from non-spin orgin and thus are due to the breaking 
of the valley or orbital degeneracy of the zero energy LL. 

As for the v — 1 QH state, the experimental data of its 
tilted field dependence, shown in the insets of Fig. 2, dis- 
plays an increasing R™™ with either a increasing in-plane 
field or a increasing normal field. While this evolution is 
consistent neither with an electron-electron or spin ori- 
gin, it is likely that the increasing shown in the inset 
of Fig. 2(a) is due to the proximity of this filling factor to 
the increasingly insulating behavior of the bilayer at the 
charge neutrality point, discussed below. However, given 
the hierarchy of the degeneracy breaking, schematically 
shown in Fig. 1(b), the mechanism underlying the v = 1 
breaking should be the same as that for v — 3, i.e. not 
from a spin origin, consistent with theoretical predica- 
tions of the formation of a QH ferromagnet in BLG at 
high magnetic fields [!, Q] ■ 

To further understand the nature of the fully lifted LL 
degeneracy, we determine the energy of the LL splittings 
by measuring i?™" at different temperatures T. Fig. 3(a- 
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FIG. 3: (a)Arrhenius plots of 7?™ n 's as a function of 1/T at 
different fields for v = 1 state, the lines with respective colors 
are the linear fits to the data points, (b) Arrhenius plot for 
v = 2 state, (c) Arrhenius plot for v — 3 state. (d)Energy 
gap AE vs magnetic field B for different filling factors v = 1 
(open triangle red), v — 2 (open circle green), v — 3 (open 
square black). The dotted lines are linear fits, whereas the 
dashed lines are square-root fits. 



c) shows log RJ xx n versus 1/T for the v = 1, 2, & 3 states. 
Since i?™ m oc exp [—(AE - 2r)/2fc B T], where AE is the 
energy gap between two subsequent LLs and T is the LL 
energy broadening. The observed Arrhenius behavior in 
these plots allows us to estimate AE at different fields 
from the slope of the line fits. Fig. 3(d) displays the field 
dependence of the activation gap for v = 1, 2, &3. Gen- 
erally, AE increases with increasing B as expected. It is 
also noted that at given B, AE V=2 > AE V = 3 , AE v —i, 
indicating that the even v states have larger energy than 
the odd v QH states in accordance with the LL symmetry 
breaking hierarchy of Fig. 1(b). As for the field depen- 
dence, we find that a \[B fit is better for v = 2 and v = 3 
states. Attempts to fit the gap evolution linearly to B 
result in a positive y-axis intercept for the v = 2 and 
v = 3 gaps and a thus a non-physical negative LL broad- 
ening r. From the y-intercept of the yB dependence, 
we can extract a physically reasonable LL broadening of 
r/fcs = 1.4 and 3.7 K for the v = 2 and 3 states, re- 
spectively [HI. Although the apparent \/S-dependence 
provides an independent confirmation of non-spin origin 
of the QH states v = 1, 2, 3, we note that the observed en- 
ergy scale AE is still too small compared to the Coulomb 
energy or even the bare Zeeman energy. For B = 15 T, 
the Coulomb energy is 370 K and the Zeeman energy 
E z = 21 K, both larger than AE V=2 = 8.6 K at this 
field. It is possible that large amounts of disorder near 
the charge neutrality point of BLG are responsible for 
such a reduced transport gap. 

We finally focus our attention on the v — QH state. 
As with SLG, the presumed v = splitting at the charge 
neutral point is not directly observable as a zero in R xx . 
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FIG. 4: Maximum resistance, measured in a two-probe con- 
stant voltage bias method, as a function of magnetic field, 
which is applied normally to the graphene plane, at a fixed 
gate voltage V S =4.1V around the charge neutrality point. The 
green dots are R^ x 's at v = using a four-probe measure- 
ment. They are taken at lower magnetic fields(_B < 20 T), 
and multiplied by a factor of 4 to match the 2 probe data. 
Upper left inset: Gatesweep of two-probe resistance at differ- 
ent fields: 19, 23, 27, and 29 T from bottom to top. All data 
are taken at 1.5K. 



Rather, R xx displays a maximum at this point, whose 
value increases with increasing B. In Fig. 4, we display 
our measurement of R xx at the charge neutrality point in 
BLG. The measured resistance shows quasi-exponential 
growth as B increases, up to ~ 10 Mil at 30 T. To avoid 
the self-heating of the graphene discussed in Ref. [ijj 
and to measure resistances > 10 Mf2 we employ a 2- 
probe AC measurement configuration with a constant 
voltage bias of 500 /zV, resulting in only ~ 10 fW of 
heating at the highest fields. The contact resistance in- 
cluded in this measurement set-up is relatively small (~ 
1 Kfi), but as an additional confirmation that it does not 
affect the behavior of the measured resistance, we cross- 
checked the 2-probe measurement using a conventional, 
current biased, 4-probe measurement at fields low enough 
that the resistance could be reliably measured with the 
voltage probes input to an amplifier with 10 MSI input 
impedance. The exponentially divergent R xx behavior at 
high magnetic field is similar to analogous measurements 
that have been performed on SLG |lll. l"l2T] , where a field- 
induced QH insulator has been proposed. We also note 
that the gate sweep (inset to Fig. 4) displays a growing 
number of local maximum in R xx , presumably due to the 
inhomogeneous distribution of these insulating states at 
high magnetic field fl9j ] . 

In conclusion, we have observed the full degeneracy lift- 
ing of the zero energy LL in bilayer graphene. Indepen- 
dent measurements of the longitudinal resistance zeros 
for the newly observed filling factors as a function of per- 



pendicular field and temperature each indicate that the 
degeneracy lifting for the v = 1, 2, & 3 splittings originate 
from electron-electron interactions. The field dependence 
of the longitudinal resistance at the v = charge neutral 
point reveals insulating behavior similar in character to 
that of single layer graphene. 
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Note added. -During the preparation of this manuscript, 
we became aware of related work with a similar conclu- 
sion from Feldman et al. 12011 . 
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